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Abstract 

We study the process e + e~ — > e + e~ $ at a high-energy e + e~ collider, where the 
missing energy arises from the radiation of Kaluza-Klein gravitons in a model 
with large extra dimensions. It is shown that at a high-energy linear collider, this 
process can not only confirm the signature of such theories but can also sometimes 
be comparable in effectiveness to the commonly discussed channel e + e~ — > 7 |5, 
especially for a large number of extra dimensions and with polarized beams. We 
also suggest some ways of distinguishing the signals of a graviton tower from 
other types of new physics signals by combining data on our suggested channel 
with those on the photon-graviton channel. 
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1 Introduction 



It is a matter of common belief in the high-energy physics community that there must be 
physics beyond the Standard Model (SM) around the TeV scale. This belief is founded in the 
fact that some of the parameters of the SM can be made stable against quantum corrections 
only by invoking new physics at the TeV scale. With this in mind, an enormous international 
effort is being poured into the construction of the 14- TeV LHC at CERN, which is expected 
to run by 2008 and come up with signals for new physics. However, as the LHC is a hadronic 
machine, it is unlikely to yield signals which are unequivocal as to the underlying theory. 
It is therefore widely deemed necessary to build a high-energy e + e~ collider to confirm the 
nature and detailed properties of this new physics. 

A high-energy e + e~ collider operating in the range 500 GeV to a few TeV must be of 
the linac type to avoid large synchrotron radiation losses. Such a machine would require a 
very high luminosity in order to be competitive with other accelerators. This is, in fact, an 
essential feature of the designpQ of proposed linear colliders. The relatively clean environment 
of an e + e~ collider allows the complete reconstruction of individual events, making it possible 
to carry out precision tests of new physics. Moreover, one has the facility of polarizing the 
beams, which can play a significant role in background elimination. It is also possible to 
run such a machine in the e~e~, ey and 77 modes. In this paper we consider the possibility 
of using a high energy e + e~ collider to detect signatures of low-scale quantum gravity in 
a model with large extra dimensions. In particular, we focus on Bhabha scattering with 
radiated gravitons. 

One of the most exciting theoretical developments of recent years has been the idea 
that there could be one or more extra spatial dimensions and the observable Universe could 
be confined to a four- dimensional hypersurface in a higher-dimensional 'bulk' spacetime[2|. 
Such ideas, which can be motivated by superstring theories, give rise to elegant solutions |S| to 
the well-known gauge hierarchy problem of high energy physics, which is just the instability 
against quantum corrections mentioned above. What is even more interesting, perhaps, is 
the suggestion [HE] that there could be observable signals of quantum gravity at current and 
future accelerator experiments, and this possibility has spawned a vast and increasing body 
of work over the past five years. This relatively new set of ideas, commonly dubbed 'Brane 
World Phenomenology', bases itself on two main principles: the concept of hidden compact 
dimensions and the string-theoretic idea of D p -branes. The simplest brane-world scenario 
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is the so-called Arkani-Hamed — Dimopoulos — Dvali (ADD) model in which there are d 
extra spatial dimensions, compactified on a d-torus of radius R c each way. Together with the 
four canonical Minkowski dimensions, this constitutes the 'bulk' spacetime. In this scenario 
the radius R c of the extra dimensions can be as large as a quarter of a millimeter However, 
the SM fields are confined to a four- dimensional slice of spacetime, with thickness not more 
than 10~ 17 cm, which is called the 'brane'. If the ADD model is embedded in a string- 
theoretic framework, the 'brane' is, in fact a Z^-brane, i.e. a 3+1 dimensional hypersurface 
on which the ends of open strings are confined[7j. However, it is not absolutely essential to 
embed the model in a string theory, and the word 'brane' or 'wall' is then used simply to 
denote the hypersurface (or thin slice) where the SM fields are confined. A crucial feature 
of this model is that gravity, which is a property of spacetime itself, is free to propagate in 
the bulk. As a result 

• Planck's constant in the bulk Ms, identified with the 'string scale', and fixed by the 
string tension a', is related to Planck's constant on the brane Mp (~ 1.2 x 10 19 GeV) 
by 

(M s ) 2+d = (4vr) d/2 T(d/2) Mp (R c )~ d (1) 

which means that for R c ~ 0.2 mm, it is possible to have Ms as low as a TeV for 
d > 2. (The normalization of reference [Sj has been adopted here). This solves the 
gauge hierarchy problem simply by bringing down the scale of new physics (i.e. strong 
gravity in this case) to about a TeV and thereby providing a natural cut-off to the SM, 
since the string scale Ms now controls graviton-induced processes on the brane. 

• There are a huge number of massive Kaluza-Klein excitations of the (bulk) graviton 
field on the brane, with masses m n = n/R c , and these collectively produce gravitational 
excitations of electroweak strength, which may be observable at current experiments 
and those planned in the near futurejH]. 

It is only fair to mention that a major drawback of the ADD model is that it creates a 
new hierarchy between the 'string scale' M s ~ 1 TeV and the size of the extra dimensions 
R^ 1 ~ 1 /ieV In fact, the huge size of the extra dimensions (compared to the Planck length) 
is not stable under quantum corrections, which tend to shrink it down until M s ~ R^ 1 ~ 
M P ~ 10 19 GeV, at which stage the original hierarchy problem is reinstated. Nevertheless, 
there are several variants of the ADD model which address this problem in various ways, 
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and some of these ideas may not be far from the truth. From a phenomenological point of 
view, it is, therefore, reasonable to postpone addressing the stability issue, and proceed to 
study the minimal ADD model and its consequences for experiment. 

The experimental consequences of ADD gravity have been mainly studied in the context 
of of (a) real graviton emission and (b) virtual graviton exchange. In the former case the 
final state gravitons in the ADD model are 'invisible', escaping the detector because of their 
feeble individual interactions (~ Mp 1 ) with matter. The final state, involving missing energy 
due to gravitons, will be built up by making an incoherent sum over the tower of graviton 
modes. In the case of virtual gravitons, the final state is built up by making a coherent sum. 
In either case, it may be shownjUIHj that, after summing, the Planck mass Mp cancels out 
of the cross-section, leaving an interaction of near-electroweak strength. 

At an e + e~ collider, the most-frequently discussedjU EI] signal for ADD gravitons 
is the process e + e~ — > 7G n , where G n is the nth Kaluza-Klein excitation of the graviton 
field. This process leads to single photon events with missing energy and momentum, and is 
expected to be among the earliest signals of ADD gravity at an e + e~ collider. In addition, 
the process e + e~ — > n + n~G n has also been studied[TT]. In this work, we have studied the 
process 

e + e~ — ► e + e~G n 

which is simply Bhabha scattering with a radiated graviton. The final state would contain 
an electron-positron pair with substantial missing energy. The reasons for taking up a study 
of this process are as follows. 

• At an e + e~ collider, final states with an e + e~ pair will be one of the first things to be 
analyzed, since this is needed for beam calibration pQ. Among such states, it should be 
a simple matter to select events with large missing energy as well. 

• Unlike e + e~ — > fi + fi~G n , this process does not suffer from large s-channel suppres- 
sion, especially if one considers final state electrons (positrons) in relatively forward 
(backward) directions. 

• Taken in conjunction with the process e + e~ — ► 7G n , this process could help to distin- 
guish signals of the ADD scenario from other physics models beyond the SM. This is 
elaborated in the subsequent discussions. 



3 



• In some ways, as we shall see, this process has signatures more distinct from other new 
physics effects than e + e~ — > jG n . 

In the following section we describe the salient features of our calculation of the signal 
process. Backgrounds and some strategies for their elimination are discussed in section 3. 
Section 4 contains the results of our numerical analysis. In section 5 we discuss how this 
process can be used to pinpoint the model if signatures of the suggested type are indeed 
observed. We summarize and conclude in the final section. 

2 The Process e + e~ — ► e + e~G n 

As gravity couples to the energy-momentum tensor, each component of the tower of ADD 
gravitons couples to all SM fields, as well as to each SM interaction vertex. The relevant 
Feynman rules can be found, for example, in Refs.jH Using these, it can be shown 
that the process e + e~ — *> e + e~G n is driven, at tree-level, by 28 Feynman diagrams. Some 
representative diagrams are shown in Fig. 1, which are obtained by 'dressing' an s-channel 
process e + e~ — > e + e~ with graviton radiation. Fig. 1(a) represents a graviton emission from 
one of the external electron legs; there will be 4 such diagrams with the graviton emitted from 
each leg in turn. Fig. 1(6) represents a graviton emission from the gauge boson propagator. 
Fig. 1(c) represents a graviton emission from one vertex; there will be another such with 
graviton emission from the other vertex. Thus Fig. 1 encompasses 14 diagrams, 7 each with 
Z and photon exchanges. In addition, there will be another set of 14 diagrams with the 7, Z 
exchange in the t-channel, making 28 diagrams in all. This latter set is absent in the case of 
e + e~ — > n + n~G n . 




Figure 1. Representative Feynman diagrams for the process e + e — > e + e G 
Considering the process 

e~(A;i,Ai) e + (k 2 ,\ 2 ) -> G n (Pi) e + (p 2 ) e~(p 3 ) 
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with a single Kaluza-Klein mode G n (pi) in the final state, we obtain the unpolarized cross- 
section 



(2) 



Ai,A2 



where the final state helicities are summed over. This cross-section is a function of the mass 
m n of the graviton excitation. 

The squared matrix element |M n (Ai, A2) | 2 for the process can be written as the sum 



|M n (A!,A 2 )p 



28 



£M«(Ai,A 2 



i=l 



2 



(3) 



where |M^(Ai, A2) | 2 represents the squared helicity amplitude arising from the ith diagram. 
In our calculation, we make use of the helicity amplitude technique to write down amplitudes 
for all these 28 diagrams and evaluate the individual terms in the above sum using the 
subroutine HELASfUJ. 

Recalling that the observed final state, namely, e + e~ plus missing energy, consists of an 
incoherent sum over the tower of graviton modes, we get 

cr(e + e~ — > e + e~+ = o~{m n ) 

n 

dm p(m) a(m) (4) 



approximating the discrete (but closely-spaced) tower of states by a continuum. The density 
of states is given byjH] 

. 2 R d c m d - 1 

P{m) = (4,r)'/* T(d/2) ' (5) 

and the integration is cut off at the kinematic limit y/s. 

As the calculation, even using HELAS, is long and cumbersome, some checks on the 

numerical results are called for. The most useful check is provided by the Ward identities 

arising from general coordinate invariance, which constitute an essential feature of any theory 

involving gravity. We can write the amplitude for the emission of any graviton in the form 

M n (X u X 2 )=T^(X u X 2 )e^>( Pl ) (6) 

where e^y(px) is the polarization tensor for the nth (massive) graviton mode. The ten- 
sor T M1/ (Ai, A 2 ) is the same for every mode, including the massless mode ej!$(pi), which is 
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the usual graviton of four-dimensional Einstein gravity. This must now satisfy the Ward 
identities 

rfT^(Ai, A 2 ) = pST^Ai, A 2 ) = , (7) 

where we note that 

28 

T^(X 1 ,X 2 ) = J2Tr^iA2) (8) 

i=l 

with % indicating the ith diagram, as above. The consistency check therefore requires a 
perfect cancellation, for each choice of Ai and A 2 , between all such terms 2 , which is highly 
sensitive to errors in signs and factors. We have also checked that if only the s-channel 
results are taken our numerical results are in close agreement with those of Ref.pi], and 
this, of course, implies that the t-channel contributions can be easily checked using crossing 
symmetry. 

In our numerical analysis of this problem, we have also considered diagrams mediated 
by gravitons in place of 7, Z. Such diagrams, in fact, are of the same order in perturbation 
theory after summation over the graviton propagators, but it turns out that the contributions 
are very small — typically less than 1% for experimentally- allowed values of Ms — due to 
the smallness of the effective graviton coupling after summation. 

It is also worth mentioning that contributions due to graviscalars in the ADD model can 
also be neglected. It can be shown |5] that diagrams with graviscalars coupling to external 
fermion legs and vertices undergo large cancellations and the final contribution is propor- 
tional to the electron mass. This leaves only the diagrams with graviscalars emitted from 
the massive Z-propagator, which is again suppressed at high energies by a factor of M|/s 
compared to the gravitensor contributions, the latter being jacked up by derivative couplings. 

3 Background Elimination 

The actual signal for the process considered in the previous section consists of a hard electron 

and a hard positron, with substantial missing energy. The SM backgrounds arise from all 

processes of the form e + e~ — > e + e~b>v, where the neutrinos can be of any flavor. At the 

lowest order, there are 96 Feynman diagrams which give rise to this final state. Some of the 

2 In fact the s and t-channel diagrams form two independent gauge-invariant sets, so that the actual 
cancellation takes place between 14 diagrams at a time. 
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principal sub-processes are listed below. 



e + e — > e ± W^v\ — > e + e v e v e 



e + e — > e + e Z — >■ e + e (z/i/) 

e + e~ — > v e v e Z — >■ z/ e i/ e (e + e~) 

e+e" -> W^+iy- -> (e + z/ e )(e-z/ e ) 

e + e~ — > ZZ — > (e + e~)(i/i?) 

-> 7*Z -> (e + e-)(z/z>) (9) 



e e 



Of course, for every process with real W/Z, there will also be many diagrams with off- 
shell particles. Taken together, all these diagrams constitute a formidable background to 
the suggested signal. A judiciously chosen set of kinematic cuts, however, enable us to 
reduce these backgrounds enormously. For example, diagrams with the neutrino pair (missing 
energy) arising from a real Z boson can be easily removed by putting a cut on the vv (missing) 
invariant mass. A similar cut on the e + e~ pairs would seem obvious but it is not advisable, 
since it would remove a large part of the signal arising from radiative return to the Z-pole 
through graviton emission. This, however, is not a serious problem, as the background 
contribution due to e + e~ — > v e v e Z —>■ v e v e {e + e~) is not very large. On the other hand, the 
background due to the process e + e~ — > e ± W T is more problematic, as most cuts tending 
to reduce this also affect the signal adversely. 

In principle, an additional source of background can be the 'two-photon' process e + e~ — > 
e + e~e + e~ where one electron-positron pair escapes detection by being emitted close to the 
beam pipe. Here one has to remember that such forward electrons, being rather hard (with 
energy on the order of 100 GeV or above), should be detectable in the end caps of the 
electromagnetic calorimeter, although their energy resolution will be rather poor£Q. Thus it 
should be possible to impose a veto on hard electrons up to within a few degrees of the beam 
pipe, whereby the above background can be virtually eliminated. Events with the forward 
electrons even closer to the beam pipef!]]] are removed via a cut on missing p?- 

Taking all these considerations into account, we focus on an e + e~ collider operating at 
a center-of-mass energy of 500 GeV (1 TeV). The initial set of kinematic cuts used in our 
analysis are as follows: 

• The final state electron (positron) should be at least 10° away from the beam pipe. 
This tames the collinear singularities arising from t-channel photon exchange. At the 
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same time, it also ensures that any background effects from beamstrahlung are mostly 
eliminated |14j. 

• The electron (positron) should have a transverse momentum p e T > 10 GeV. 

• We demand a missing transverse momentum p™ lss > 15(25) GeV for y/s = 500 GeV 
(1 TeV). As mentioned above, this also helps in reducing the two-photon background. 

• The tracks due to electron and positron must be well-separated, with AR > 0.2, where 
R = y/Ar] 2 + A<p 2 , in terms of the pseudorapidity rj and the azimuthal angle (p. 

• The opening angle between the electron and positron tracks is required to be limited 
by 5° < 6 e + e - < 175°. This ensures elimination of possible cosmic ray backgrounds, 
and also ensures sufficient missing energy. 

• The missing invariant mass M miss should satisfy the cut |M miss — Mz\ > 10 GeV. This 
important cut eliminates many SM backgrounds in which two final state neutrinos 
arise from a Z-decay. 
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Figure 2: Distribution in missing invariant mass M m - lss for s/s = (a) 500 GeV and (b) 1 TeV. In addition 
to total signal and background cross-sections, we also show the signal due to s- channel diagrams only, which 
is identical with the cross-section for e + e~ — > fj, + fj,~G n . 

Incorporating all these cuts into a Monte Carlo event generator, we obtain the distribution 
in missing invariant mass shown in Figure 2. It is clear from the figure that the missing 
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invariant mass of the signal is much harder than that of the background. This is easy to 
understand if we realize that the missing mass is a measure of the energy of the emitted 
graviton. Now it is well-known that the higher the energy, the more strongly the graviton 
is coupled to matter and consequently the higher the production cross-section. Moreover, a 
higher energy leads to a higher density of states, and this further enhances the cross-section. 
In Figure 2, we have also shown the value of the signal when only s-channel diagrams are 
considered, which would be the case if we were to study processes like e + e~ — > [i + [i~G n or 
e + e~ — > T + r~G n . Comparison with the signal of interest makes it obvious that there is a 
substantial contribution to the Bhabha scattering signal from the i-channel, which makes 
it more viable as a signal of ADD gravity than the other two. In fact, the emission of a 
high-energy graviton (corresponding to large missing invariant mass) reduces the energy of 
the underlying Bhabha scattering process, causing a large enhancement due to the Z-pole. 
This is the reason why the signal for large M miss is completely dominated by the s-channel 
contribution. At the same time, when the graviton is relatively soft, the s-channel diagrams 
are strongly suppressed, but the t-channel contributions are not; as a result the signal is 
dominated by the t-channel contribution for low M miss . We can therefore obtain a clear 
separation of signal from background by imposing a cut 

M miss > 350 (450) GeV for = 500 GeV (1 TeV) 

These choices are optimal, after taking into account the fact that the signal falls rapidly with 
increasing M s and d. 

An important strategy for background reduction is the use of the beam polarization 
facility at a linear collider. If the electron (positron) beam has a right (left) polarization 
efficiency V e (Vp), the cross-section formula corresponding to equation El is obtained by the 
replacement 

£ \M n (X 1 ,X 2 )\ 2 — (l + V e )(l-V p )\M n (+,+)\ 2 + (l + V e )(l+V p )\M n (+,-)\ 2 

Ai,A2 

+ (1 - V e ){l - V p )\M n {-, +)| 2 + (1 - P e )(l + V p )\M n (- -)\ 2 

(10) 

Typical values for the polarization efficiencies e T used in this analysis are V e = 0.8 and 
V p = 0.6, which tend to favor the first two terms on the right side of the above expression. 
The use of polarized beams leads to a drastic reduction in the SM background, chiefly because 
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it causes suppression of the VF-induced diagrams. The numerical effects are presented in the 
next section. 



4 Results and Discussions 

After numerical evaluation of signal and background, we find it convenient to present our 
results, using the significance S/yB, where S(B) = C<ts(b)-, for given values of d and Ms, 
Here <Js(B) denotes the cross-section for the signal (background), while C is the integrated 
luminosity. Our numerical results are presented assuming £ = 500 fb _1 ; however, it is a 
simple matter to scale the significance values for different values of luminosity. 

In Figure 3(a) we show curves showing the variation of the significance with the string 
scale Ms, for d = 2, 3, 4, 5 and 6 respectively, for = 500 GeV. Figure 3(6) shows a similar 
plot for yfs = 1 TeV. In both cases, we have considered unpolarized electron and positron 
beams. 
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Fi gure 3: Illustrating the variation of significance S/yB with Ms, considering unpolarized beams for (a) 
v/i = 500 GeV, and (b) y/s = 1 TeV. 



As has been already mentioned, the use of beam polarization improves changes of detect- 
ing the signal quite dramatically. This is illustrated in Figure 4, which is similar to Figure 
3, except that the polarization efficiencies have been taken to be V e = 0.8 and V p = 0.6 for 
the electron and positron respectively. 
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Fi gure 4: Illustrating the variation of significance S/y/B with Ms, considering polarized beams for (a) 
Vs = 500 GeV, and (b) y/s = 1 TeV. We take V e = 0.8 and V p = 0.6. 



From the above graphs, it is straightforward to find out the maximum value of M s that can 
be probed at the linear collider for any given value of d, both for polarized and unpolarized 
beams. In Table 1, we present the values of M s that one can probe at 99% confidence level. 





d = 


2 


3 


4 


5 


6 


0.5 


polarized 


4.53 


2.79 


2.01 


1.60 


1.35 




unpolarized 


3.16 


2.09 


1.59 


1.30 


1.12 


1.0 


polarized 


7.40 


4.63 


3.41 


2.76 


2.36 




unpolarized 


5.36 


3.57 


2.74 


2.28 


2.36 



Table 1. 99% C.L. discovery limits on the string scale Ms for an ADD scenario using (graviton) radiative 
Bhabha scattering, for different numbers of extra dimensions. An integrated luminosity of 500 /6 _1 has been 
assumed. All numbers are in TeV. 

If one compares these with the corresponding reach of the processjH EH] e + e~ — ► 
7+ ft, one will notice that search limits are approximately of the same order. However, 
the comparison should not be made too literally, because of several reasons. First, the 
polarization efficiencies in our case, while matching with those of [TP] , are different from 
those of |5 , while the center-of-mass energy at which reference ^U] has calculated the effects 
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is slightly different from ours. In addition, one needs to match the event selection criteria 
more carefully for a full-fledged comparison. Finally, it should be borne in mind that in 
our analysis we have adopted a normalization of the string scale M s , which, though widely 
used, is not necessarily uniform in the literature. However, in spite of such non-uniformities, 
Figures 3 and 4, together with Table 1, can perhaps be taken as faithful indications of the 
fact that the predictions on radiative Bhabha scattering are comparable to those on the 
photon-graviton channel, so far as the limits of the probe on the string scale in ADD models 
are concerned. 

5 Pinning Down the Model 

A few comments are in order on how to distinguish graviton signals involving missing energy 
from similar ones arising from other kinds of new physics. For example, let us consider the 
well-known process e + e~ — > ^yG n , leading to a single-photon-plus-missing-energy signal. Such 
signals can be obtained in several other models. A model with an extra Z' boson could lead 
to a process e + e _ — > jZ' with the Z' decaying to neutrinos. Similarly, in supersymmetric 
models, a process like e + e~ — > 7X1X1 (where x? is the lightest neutralino) would also lead 
to a single photon signal. Although the Z 1 signal can be differentiated by looking for a 
resonant peak in the missing invariant mass, the supersymmetric signal is more difficult to 
disentangle. 

The situation for the process considered in this paper is somewhat more encouraging. It 
is true that both types of new physics considered in the preceding paragraph can produce 
the same signal. Processes of the form e + e~ — > e + e~Z', with the Z 1 decaying to neutrinos 
are one possibility. In supersymmetry, too, it is possible to pair-produce either sleptons or 
charginos, followed by decays to electrons (positrons) and invisible particles. The Z' signal 
can again be differentiated by a peak in the missing invariant mass. The supersymmetric 
signal, in this case, always arises from production of a pair of real sparticles, which should 
emerge in opposite hemispheres and, if light enough, will be highly boosted. As a result, the 
observed electron and positron should lie, most of the time, in opposite hemispheres. For the 
graviton signal, however, there is no such correlation; in fact the distribution in the opening 
angle between e + and e" turns out to be practically flat. This distinction clearly does not 
work very well when the produced sparticles are heavy. 

A more elegant way of distinguishing these signals from those other forms of new physics 
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<*e + e -> e + e G„[P b ] <*e + e -> e + e G„[P b ] 

Fi gure 5: Correlation plot showing the polarized cross-sections for the processes e + e — > jG n and e + e — > 
e + e~G n for different values of d. Each curve is generated by varying Ms for (a) y/s = 500 GeV, and (b) 
yfs = 1 TeV. Shading indicates the region where the theoretical calculation is unreliable. 

is to compare the cross-sections arising from both the processes e + e~ — > 7G n and e + e~ — > 
e + e~G n , both of which are determined by the two parameters d and Ms, and hence, will have 
some correlation. In Figure 5, we have plotted the cross-sections a{^G n ) versus a{e + e~G n ) 
for different values of d. In view of the discussion in the previous section, we would like to 
note that here the rates for e + e~ — > jG n have been computed with the same normalization 
of the string scale as that used for the Bhabha scattering process. Each curve corresponds 
to variation of Ms over the range allowed by current experimental constraints. Once exper- 
imental data are available from a high energy e + e~ machine, we should be able to pinpoint 
a small region (ideally a point) on the graph(s) in Figure 5. The position of this would 
immediately direct us to the value of d; comparison of the cross-sections would now yield 
a measurement of the value of Ms- The identification of the model parameters should be 
unambiguous since the curves for different values of d do not intersect except at the origin. 
In any case, the area in the vicinity of the origin is limited by the search limits indicated 
in Table 1 since it corresponds to very high values of Ms- The shaded region corresponds 
to Ms = y/s beyond which the theoretical calculations are unreliable. It should, of course 
be remembered that further information on the fundamental parameters of theory can be 
extracted from various kinematic distributions which are sensitive to these parameters. 
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6 Summary and Conclusions 



In the paper we have considered the process e + e~ — > e + e~G n , which is essentially Bhabha 
scattering with radiated gravitons in the ADD model. After identifying suitable event selec- 
tion criteria, we find that this process can act as an effective probe of large extra dimensions 
at a high-energy e + e _ collider, especially with a center-of-mass energy of order 1 TeV and 
with polarized beams. The string scale Ms that can be probed in this channel is found to 
be comparable to that which is accessible through the alternative process e + e" — ► 7G n . It is 
also shown that a study of the (graviton) radiative Bhabha scattering process may provide 
some further handle on the essential characteristics of ADD-like theories. And finally, taking 
a cue from the wisdom that, in coming to any conclusion on new physics possibilities, it is 
always advantageous to have more than one type of data, we have demonstrated how our 
predictions can be combined with those on the photon-graviton channel to obtain rather 
trustworthy revelations on models with large extra dimensions. 
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